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Abstract 

We study the low temperature dynamics of a single flux line in a bulk type-II 

superconductor, driven by a surface current, both near and above the onset 

of an instability which sets in at a critical driving. We found that above the 

critical driving, the velocity profile of the flux line develops a discontinuity. 
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Dynamics of a driven elastic string have attracted much recent attention [0,0. While 
most of the work has been focused on the interesting physics of pinning-depinning transitions 
in the case of bulk driving, the paper by Tang, Feng, and Golubovic [0] studied the case 
of a surface-current-driven flux line in a bulk type-II superconductor. They found a novel 
instability of the flux line motion at large driving currents. The instability sets in at a critical 
driving, where the line loses its steady state motion and (presumably) will be stretched 
longer and longer. Their finding depends crucially on the boundary condition they use. 
Physically, the surface driving current is within a boundary layer of thickness A, where A is 
the penetration depth. The boundary condition used in Ref. || is somewhat equivalent to 
taking the limit A — > in a plausible but uncontrolled way. Since the instability sets in at or 
near the boundaries, it is necessary to examine the situation carefully using a more physical 
boundary layer. Also, it is important to see what happens when the driving current is larger 
than the critical driving - a question which can not be addressed by using the boundary 
condition in Ref. 0. 

In this paper we analyze the flux motion with the more physical boundary layer Lorentz 
driving force. We first use the method of matching asymptotic expansions to study the steady 
state solutions. The lowest order matching condition justifies the form of the boundary 
conditions used in |2| and gives the relation of the driving force to the current. We then 
study, both numerically and analytically, the complete equation below and above the onset 
of instability. 

Let us first derive the equation for the flux line motion which involves the Lorentz force 
as a term in the equation, as opposed to just a boundary condition. As we will be mostly 
interested in fairly large driving forces, we neglect pinning effects. The Lorentz force on a 
flux line is just F = - J j x hdsdA where s is the arclength along the flux line and dA a 
section of infinitesimal area transverse to the flux line. If the applied current, j, is slowly 
varying in the direction transverse to the line, then the integration in these coordinates may 
be carried out to give 
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F = ^/jx ids, (1) 

where 0o is the flux quantum, and t is the unit tangent vector in the direction of the local 
magnetic field (arclength is taken to be increasing in the direction of the magnetic field). The 
exact form of the current depends on the geometry of the sample, however it is known that 
the magnitude of the applied current drops exponentially with distance from the boundary 
of the sample. 

For simplicity we model the dynamics of a single flux line as a two dimensional problem, 
defined by its shape function r(s,t), or where the parameterization is well defined, y(x,t) 
(see Fig. |I|). The applied field is in the negative x-direction, and the applied current in the 
negative z-direction, thus giving the driving Lorentz force predominantly in the y-direction. 
The Lorentz force per unit length is then 

dF 0o . f (x-L/2\ (-x-L/2\\ 

where jo is the current density at the surface, the sample boundary is at x = ±L/2, and 
n = z x t is the local unit normal vector of the flux line. 
The equation of motion for the flux line becomes 
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The left-hand side, the viscous damping, and the first term on the right, the normal force 
due to line tension, are the same as in |J, and the last term is the Lorentz force, from 
Eq. (Q). 7 is the damping coefficient (7 ~ 27r gf°2 p from the Bardeen-Stephen model ||, with 
£ the coherence length, and p n the normal state resistivity), and a is the line tension, given 
approximately by -^-47r^ 2 ln(«;), with H c the critical field and n = A/£ the Ginzburg-Landau 
parameter. K is the curvature, and we have the relations t = d s r and d s t = Kn. 

In cases where the tangent vector never becomes vertical (ie. dy/dx remains finite), 
Eq. (Q) can be rewritten in terms of the x and y coordinates of r, now reparameterized by x. 
Note that a displacement of An in the direction of the normal is related to a displacement 
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Ay = Anyl + (dy/dx) 2 and the curvature in terms of x and y is given by K = [1 + 
(dy / dx) 2 }~ 3 ^ 2 (d 2 y / dx 2 ) . Thus we get the equation 
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We now examine the steady state solutions of Eq. (^). Steady state implies v = dy/dt is 
constant, which allows us to rewrite (HI) as a first order equation in the sine of the tangent 



angle, 6. Setting / = and w = sin8 = jf-fyl + (dy/dx) 2 we get 
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For x far from the boundaries (i.e. \x ± L/2\ ^> A) the driving term is negligible and the 
equation becomes 



dx 



~v)y/l- W 2 , 
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which has solution 
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Now we examine the solution near the boundary at x = L/2. The coordinate appropriate 
in this region is r) = — x ~^ 2 - In terms of 77, our Eq. ([]) becomes 



dwj 
dr] 



7 



a 



v \Jl- w 2 + fe-^, 



(8) 



If we expand Wi in powers of (^f) A as Wi = + (^.v) Xw^ + ... we obtain a series of 
equations for the w^. The first two of these equations are 
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Assuming an applied field perpendicular to the boundary, these have the solution 
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If we expand w Q about x = L/2 (77 = 0) and iUj for large rj we get 
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I-/ 2 17. 
Matching w Q to U7j gives, to order 
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Note that the velocity found in || is obtained by dropping the term of order ^ on the right 
hand side of (|l^) and thus is the zeroth order of our asymptotic solution. This matching 
procedure is illustrated in Fig. 0b which shows w = sin6 as a function of x for / = 0.9. 
The solid line is a steady state numerical solution, and the broken and dashed lines show 
the inner and outer solutions, respectively. We see that the outer and inner solutions agree 
very well with the numerical result within their respective domains of validity. A composite 
solution, valid on the whole domain, can be formed by adding w Q and Wi and subtracting 
their common part from Eq. (|XT|) . This is indistinguishable from the numerical solution in 
Fig. |b. 

The numerical solutions shown in Fig. ^| were produced from solutions of Eq. @. This 
was chosen, rather than Eq. (|j) in x-y coordinates, due to problems arising in the continuity 
of dy/dx and the diverging values of dy/dx found at large values of / (see below). As our 
equation involves the position vector, r, explicitly, we must evolve a set of vectors {r(s)} of 
positions along the curve (as opposed to, for instance, following the curvature). We solve 



Eq. (0) using a finite-difference approach. The viscous term, jt, and the curvature term, 
Kn = r ss , can be dealt with using a Crank-Nicholson type approach for diffusive equations. 
This yields two, x and y, tridiagonal systems linked only at the boundaries. The Lorentz 
force in (|J) is then dealt with in an semi-implicit manner. The system is remeshed at each 
time step to preserve point spacing in regions of high curvature. 

For a specific case, we take a sample width, L, of 100A and measure the velocity in 
the unit of er/7. Fig. |]a shows the line shapes for / = 0.2 to / = 1.1. We see that the 
slope remains fairly small within a penetration depth, A of the boundary, consistent with 
the assumptions for Eq. ([!]). Also, the analytic solution of Ref. |§ starts to deviate from 
our numerical solution near the boundary for large /. 

Fig. |3] shows v as a function of /. The crosses are from steady state numerical solutions, 
the dashed line is the zero'th order matching condition from 0, and the dotted line (for 
the region / < 1) is from Eq. ([12]). The zero'th order solution suggests that as / — > 1, 
v — ► v max = na/'-fL, implying 9 — > 7r/2; i.e. the flux line "wets" the boundary. The more 
accurate expression, Eq. (|T2"D, suggests that v — > arcsin(l — (^)(1 + 7r/2)) < v max so 
the flux line does not become vertical as / — > 1 (see also the numerical solution in Fig. [|). 
What then does happen for / greater than one? As we shall see below, the flux line becomes 
vertical (9 — > it/ 2) at an interior point, but not until / = /* = 1.07623 for our sample case 
where X/L = 0.01. Above /* the speed of the flux line develops a discontinuity, becoming 
piecewise constant. 

Note that in the above analysis wf^ > 1 for f > 1, so it can not be extended to the region 
where / > 1. This problem can be remedied by adjusting the arbitrary constant in so 
that wf^ does not exceed one. This means that will no longer satisfy the boundary 
condition w^l^o = 0. We can, however, adjust the constant in to compensate for this 
discrepancy so that wf 1 + (-^ L )w^ = at 77 = 0. This results in a solution to Eq. (|), for 
/ > 1, of 
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Expanding this Wi for large 77 gives 
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Matching this to w a at x = L/2 gives, to order ^ 
2a 
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Note that for / = 1, v = j?- arcsin(/ — (^)(1 + 7r/2)), the same result as taking / — > 1 in 
Eq. (pD. Eq. (||) has only real solutions for f < f* = 1.07623 (for L = 100A). It suggests 
that the instability should occur at / = /* where v = v max = ira/'yL. Eq. fllS]) is shown as 
the continuation of the dotted line for 1 < / < /* in Fig. ^|. For general X/L, f* is found 
as the root of Eq. (|I~5| ) for v = v max = na/^L. For small X/L, the case we are interested in, 
this root is 

f* « l + y(7r + 7r 2 /4) (16) 

We see that as X/L — >• 0, /* 1. 

The question now arises as to what happens above /*. Fig. ^a shows the numerical 
evolution of the flux line shape for / = 1.1, just above the transition, and / = 1.5. There 
are two important things to note in this figure. First, the flux line is approaching a vertical 
asymptote at about x = 40.5 in what seems to be an asymptotic manner (i.e. the flux line 
does not become vertical in a finite amount of time). Secondly, the portion of the flux line 
to the boundary side of this vertical asymptote has a constant shape, implying it is moving 
with a constant speed in the y- direction. This last observation can be verified by applying a 
finite difference approximation to Eq. (f|) to compute dy/dt for the points on the flux lines 
of Fig. [|a. The result of this computation is shown in Fig. |]b. We see from this velocity 



profile that, indeed the speed is constant in the boundary layer, but that a discontinuity has 
developed in the velocity profile! The constant speed of the boundary layer can be deduced 
as follows. 



This is shown as the dotted line in Fig. ||] for / > /*. Comparison of this speed to the speed 
obtained in the numerical simulations shows excellent agreement. 

Now, what about the outer solution, w Q l We see from Fig. |] that as time progresses 
the speed of the outer solution approaches a constant value and that the position of the 
discontinuity in the velocity profile (or the vertical asymptote in Fig. ^a) seems to approach 
a fixed value. The location of the vertical asymptote and the asymptotic (large time) speed 
of the inner solution are quite related. Requiring that w Q — > 1 at the vertical asymptote 
gives the speed of the outer solution. 

In conclusion, we have studied the flux line motion, in particular the dynamical instability 
found in Ref. [Q], using a more physical boundary layer driving. The boundary condition 
used in [[| is consistent with our zeroth order (in A) asymptotic matching. The analytic 
solution of Ref. is quantitatively valid for / < 0.8. For larger /, the deviations both in 
line shape near the boundary and in the velocity are significant. We have shown that the 
instability occurs at / = /* = 1 + (\/L)(n + 7r 2 /4) where the flux line starts to loose a 
steady state motion. We have observed numerically that above this instability the flux line 
velocity profile develops a discontinuity. This instability has a clear mark on the "I — V" 
curve, Fig. [| that is a sharp upward turn at /*. As pointed out in Ref. [Q], this instability 
should also occur in dense flux line systems. 



Eq. (|TJD gives the constant asymptotic value of Wi for large r\. Above /* this asymptotic 
value can only be one (wi = sin 6 < 1), as the flux line becomes vertical. So setting (|TJj) to 



one gives the speed of the inner solution, v^, for / > /* as 




(17) 
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FIGURES 

FIG. 1. Sketch of the cross section of the sample. The driving current is near the sample 
surfaces. 

FIG. 2. (a) Steady state flux line profiles for, from top to bottom, 7=0.2, 0.4, 0.6, 0.8, 1.0, and 
1.1. (b) Matched asymptotic expansions for / = 0.9: numerical solution (solid line), w (dashed 
line), and u>i (dot-dashed line). 

FIG. 3. The velocity of the flux line as a function of driving force: numerical simulations 
(crosses), Eqs. (12), (15), and (16) (dotted line), and Eq. (10) of Ref. [1] which is also the zeroth 
order of Eq. (12) (dashed line). 

FIG. 4. (a) The shape of the flux line at different times for / = 1.1 (solid lines) and / = 1.5 
(dotted lines). The solutions which extend further down the plot are at later times, (b) Velocity 
profiles of the flux lines of (a) at the latest times shown. 
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